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lONIZATlON  OF  A  HIGH  ENERGY  NEUTTIAL  BEAM 

Propagating  in  the  Ionosphere 


Shu  T.  Lai 

Phillips  LaboratoryAVSSI,  Hanscom  AFB,  MA  01731 


Abstract 

By  using  modeled  stripping  ionization  cross  sections  and  standard 
atmospheric  density  models,  we  calculate  stripping  ionization  probabih'ty  of 
high  energy  (MeV)  neutral  beams  propagating  in  the  earth’s  atmosphere. 
Numerical  results  as  functions  of  altitude,  beam  angle  and  beam  energy  are 
presented.  We  also  consider  Alfv6n’s  critical  ionization  velocity  (CIV)  process 
as  an  alternate  ionization  pathway.  If  CTV  occurs  as  the  neutral  beam 
propagates  through  the  geoplasma  across  the  ambient  magnetic  field,  it  may 
be  an  important  pathway  to  rapid  ionization  halting  the  beam.  We  conclude 
that  the  parametric  conditions  are  unfavorable  for  CIV  to  occur. 

I  Introduction 

When  a  high  energy  (MeV)  neutral  beam  propagates  in  the  earth’s  ionosphere,  the 
beam  particles  interact  with  the  particles  and  fields  in  the  ionosphere.  If  the  beam  becomes 
ionized  as  a  result  of  the  interaction,  the  beam  may  be  halted  since  the  product  beam  ions 
gyrate  around  the  ambient  magnetic  field  lines.  The  cross-sections  of  electron  impact 
ionization  and  charge  transfer  are  insignificant  at  Mev  energies.  We  do  not  consider  nuclear 
reactions  because  they  do  not  occur  at  energies  below  30  MeV.  Scattering  dispersion  of  the 
beam  can  degrade  beam  energy  but  is  unrelated  to  ionization.  The  dominant  ionization 
process  for  a  MeV  neutral  beam  propagating  in  the  ionosphere  is  probably  stripping: 


N  +  A-N*  -*-A  +  e' 


(1) 


where  iV  is  a  generic  beam  particle  and  M  a  generic  atmospheric  neutral  particle  or  ion. 
Neutral  particles  are  orders  of  magnitude  more  abundant  than  ions  in  the  ionosphere. 
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II.  Theoretical  Model 


Measurements  of  the  cross  sections  a  of  stripping  of  hydrogen  H  incident  on  various 
gases  have  been  reported  [1,2,3, 4].  An  empirical  formula  a(E)  of  H  is  given  in  ReL[5]: 


o(£)-o^ 


(2a)Q(£-/)* 


(2) 


where  0^=10'^^  cm^,  n=0.75,  E  is  the  energy  (in  keV)  of  the  beam,  Z,  v,  J,  and  a  are 
parameters  of  the  ambient  species,  and  I  the  beam  ionization  energy  (in  keV). 

The  density  /(z)  of  a  neutral  beam  of  energy  E  propagating  from  altitude  z^  to 
altitude  z  is  modeled  simply  as 

I(E,z)  -  /(£.z,)exp{-/^*53.  [n.(z)o, (£)]dj/cos0}  (3) 


where  6  is  the  angle  between  the  beam  and  the  vertical  line  (8=180°  when  the  beam 
propagates  downwards).  We  take  the  summation,  i=l  to  5,  over  the  5  most  abundant 
atmospheric  species  O,  O2,  N2,  He,  and  H.  The  survivial  probability  F(£iz)  of  the  beam 
undergoing  stripping  is  given  by 

P(E.z)-l--^  (4) 

For  atmospheric  densities  we  use  the  Stein- Walker  model  [6]  for  simplici^. 


III.  Results 

The  survival  probabilities  of  two  downward  propagating  neutral  H  beams  (2  and  250 
MeV)  are  calculated  (Fig.l).  The  higher  energy  beam  can  survive  longer  until  it  reaches 
about  110  km.  where  it  is  completely  ionized.  The  survival  probability  of  a  2  MeV  beam 
propagating  horizontally  at  various  altitudes  through  various  distances  (1km  to  1000  km)  is 
also  shown  (Fig.2). 
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Figure  1.  Survival  probabilities  of  2  MeV  and  250  MeV  neutral  hydrogen 
beams  propagating  vertically  downwards  in  the  ionosphere. 


Figure  2.  Survival  probability  of  a  2  MeV  neutral  hydrogen  beam 
propagating  horizontally  in  the  ionosphere. 
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IV.  CIV  Criteria 

Alfven’s  [7]  critical  ionization  velocity  (CIV)  suggests  that  when  a  neutral  gas  and  a 
magnetized  plasma  travel  relative  to  each  other  with  a  velocity  exceeding  a  critical  velocity 
=  V(2e<f>/M),  rapid  ionization  occurs.  M  is  the  mass  of  a  neutral  particle  and  e0  the 
ionization  energy.  When  an  ion  beam  travels  across  the  ambient  magnetic  field,  beam- 
plasma  interaction  occurs.  As  a  result,  the  electrons  form  a  plateau  tail  distribution.  Some 
electrons  in  the  tail  may  be  energetic  enough  to  ionize.  For  a  review  of  CIV,  see,  for 
example,  Ref.  [8]. 

Carini  et  al  [9]  questioned  whether  CIV  can  occur  in  MeV  neutral  beams.  If  CFV 

occurs,  it  could  be  a  rapid  path  to  ionization  halting  the  neutral  beam. 

•  » 

For  CrV  to  occur,  it  is  necessary  [10]  that  an  electron  has  to  ionize  at  least  once 
before  it  leaves  the  interaction  region.  This  criterion  requires  tv  >  1  where  t  is  the 
electron  transit  time  and  v  the  ionization  frequency.  For  a  narrow  low  density  beam,  both 
T  and  V  are  small. 

It  is  also  necessary  that  the  contact  time  of  a  neutral  beam  pulse  with  an  ambient 
magnetic  field  line  should  be  longer  than  the  electron  energization  time  [11]. 

t,>r,  (5) 

For  a  MeV  H  beam  pulse  of  length  Im,  -O.ljtxs.  Taking  -30/ [12]  where 
is  the  lower  hybrid  frequency,  one  finds  that  in  the  ionosphere  -10ms,  which  exceeds  Tj^. 
Electromagnetic  modes  may  be  excited  when  the  beam  velocity  V  exceeds  [13]. 

+  (6) 

where  jS  ‘•nkTIB^  and  is  the  Alfven  velocity.  These  modes  drain  energy  and  may 
suppress  CIV.  Thus,  for  CIV  to  occur,  the  beam  velocity  criterion  is 

V^>V>V^  (7) 


B 
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N 


Since  -  10*^  and  (O)  -  7xl0^km/s  in  the  ionosphere,  it  appears  that  CIV  is  unlikely 
because  the  beam  velocity  V  (-10^  km/s)  exceeds  (-7x10^  km/s). 
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